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Abstract 

Transport theory for an interacting fermionic system is reviewed and 
applied to the chiral Lagrangian of the Nambu-Jona-Lasinio model. 
Two expansions must be applied: an expansion in the inverse number 
of colors, l/Nc, due to the nature of the strong coupling theory, and 
a semiclassical expansion, in powers of h. The quasiparticle approxi- 
mation is implemented at an early stage, and spin effects are omitted. 
The self-energy is evaluated, self-consistently only in the Hartree ap- 
proximation, and semi-perturbatively in the collision integral. In the 
Hartree approximation, 0((1/A'^c)'^)! the Vlasov equation is recovered 
to 0{h^), together with an on- mass shell constraint equation, that is 
automatically fulfilled by the quasiparticle ansatz. The expressions for 
the self-energy to order 0{{1/Nc)) lead to the collision term. Here 
one sees explicitly that particle-antiparticle creation and annihilation 
processes are suppressed that would otherwise be present, should an 
off-shell energy spectral function be admitted. A clear identification 
of the s, t and u channel scattering processes in connection with the 
self-energy graphs is made and the origin of the mixed terms is made 
evident. Finally, after ordering according to powers in h, a Boltzmann- 
like form for the collision integral is obtained. 



*Dedicated to Professor R.H. Lemmer on the occasion of his 65th birthday 
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I. INTRODUCTION 



The description of physical systems that are out of thermal and chemical equi- 
librium has its origins in the early work of Boltzmann in 1872 which describes the 
development of a system via an equation for the single particle distribution func- 
tion. The formal derivation of non-relativistic transport equations from a specific 
field theory can be traced back to the work of Kadanoff and Baym [Q], but which can 
be succinctly stated in a closed time path (CTP) formalism that is independently 
due to Schwinger and Keldysh P| To date, this formalism has been applied to 
many non-relativistic problems, but also to some in its relativistic generalization. 
There are several excellent books and reviews available which elucidate the 

techniques, both from an operator or path integral point of view, to which we refer 
the reader. 

Recently there has been an additional surge of interest in non-equilibrium for- 
mulations of field theory from the nuclear physics community, as a need exists for 
describing heavy-ion collisions within such a framework. A seminal article on this 
topic is that of Mrowczyhski and Heinz p|, who, independently of Schonhofen et al. 
, have derived transport equations for the Walecka model. The derivation of such 
transport equations is in fact more general than appears on first sight, and the basic 
equations of motion for the Green functions, as they occur in Ref. are applicable 
to any field theory. The major difficulty in the fermionic sector lies in the fact that 
(a) there are in general 16 x 2 equations to be solved, and (b) should one in any 
event introduce a simplifying assumption such as the quasiparticle approximation, 
one has to perform a complicated evaluation of the self-energy that goes beyond 
the mean field approximation, if one wishes to describe interesting physics. While 
this has often been done (see for example, Ref. @,|l3])5 it is often unclear how the 
final collision integral is related to scattering diagrams in a strict Feynman sense, 
although one knows heuristically that the Boltzmann equation must emerge from a 
semi-classical approximation. 

The purpose of this paper is twofold. On the one hand, we wish to make clear 



Some work that already introduced the necessary concepts by Bakshi and Mahanthappa 
was already published in Ref. Q 
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how the coUision term in a specific model is related to the Boltzmann equation 
via Feynman (and not heuristic) diagrams in evaluating the scattering cross section 
that appears in this expression, and also to explicitly demonstrate the limits of the 
quasiparticle assumption in that it supresses vacuum fluctuations and particle pro- 
duction and annihilation processes that would be present if the on-shell condition 
were relaxed. Although the calculation is performed for a specific model, its gen- 
erality is evident. A transport calculation for QCD, or any other theory involving 
relativistic fermions, for that matter, should also display the structure given in this 
paperQ Thus this calculation gives a working guide for the general construction of 
the collision term. 

The second purpose of this manuscript is to develop the transport theory for a 
fermionic chiral Lagrangian, in our case, the minimal version of the Nambu-Jona- 
Lasinio (NJL) model |J2|,|13|. This specific choice is motivated by the following 
physical considerations. Lattice simulations of QCD, which are equilibrium studies 
of this field theory, show a chiral phase transition at a fixed value of the temperature 



14], a feature which can be simply modelled by the NJL Lagrangian, see for example 



T5[] . Unfortunately, due to the difficulties inherent in QCD itself, there appears to be 



no means of observing such a phase transition in an equilibrium system and extract, 
for example, critical exponents. Rather, the experimental search for clues to a phase 
transition that has been conceived, is to proceed via heavy ion collisions, in which 
both high temperatures and densities are to be reached. Much effort has been 
expended in this direction and will continue to be so during the early course of the 
next century, particularly at RHIC and CERN. In addition, the simple equilibrium 
picture of the static lattice simulations is inadequate for a complete description of 
such heavy ion experiments. We note that, from the theoretical side, there is at 
present no single conclusive observable that might indicate the occurrence of the 
chiral phase transition, given that one is dealing with a non-equilibrium system^. 



This is not evident, for example, in the diagrams shown by Geiger [11| who omits 
exchange graphs. 

^ Some research has been undertaken on the basis of static calculations that imply that 
substantial changes my occur in the mesonic spectrum, in particular with regard to the p 
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Thus we may speculate as to whether a concrete signal of this phase transition 
can be obtained from a non-equilibrium formulation of the problem or transport 
theory, for a simplified field theory that is based on an underlying chiral Lagrangian 
and which therefor displays the same symmetry structure as QCD. Here one may 



surmise whether phenomena such as critical scattering i.e. divergence of the 
scattering cross-sections at the phase transition, may leave an obvious signal, or 
simply the fact that the dynamically generated quark mass is a function of space 
and time may leave some visible remnant of a phase transition. 

Our starting point is thus the Nambu-Jona-Lasinio (NJL) model, that intrinsi- 
cally contains only quark degrees of freedom that interact via an effective two body 
(in SUf{2)) and, in addition, three body interactions for three flavors. Mesonic 
states are constructed as collective excitations in this model. No explicit gluonic 
degrees of freedom are included. For our purposes, it is sufficient to study the sim- 
plest form of this model that displays the features of chiral symmetry, and we will 
thus restrict ourselves to the simplest form of Lagrange density that does so. 

For the derivation of transport equations, the Schwinger-Keldysh formalism 
for non-equilibrium Green functions is applied, in what is by now a standard ap- 
proach We follow closely the approach of Ref. [|]. The associated equations 
of motion for the Wigner functions are then written down in a general manner that 
pertains to any fermionic theory. The difficulty in solution, on the other hand, as 
has been pointed out, has its origins primarily in the complication caused by the 
fact that one is dealing with relativistic spinors |T^ . 



In this paper, we construct the transport and constraint equations for the NJL 
model Lagrange density that we have chosen. There are four esssential ingredients 
to this calculation. 

(i) We choose to make a quasiparticle expansion of the Wigner function at a 
rather early stage of the derivation, and (ii) we do not allow for spin effects. This 
reduces the independent functions to simply be the distribution functions for quarks 
and antiquarks for which equations of motion are derived, whose ingredients are the 
self-energies. These two assumptions, while not necessary in principle, reduce the 



mass. Such changes, however, are inferred from model calculations, usually together with 



additional assumptions |14], and they differ according to the approach taken. 
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complexity of the manipulations tremendously. They also correspond to what is 
at present possible in numerical calculations, (iii) Noting that the NJL model is 
a strong interaction theory, one requires a non-perturbative expansion of the self- 
energies in the inverse number of colors l/N^. |T8| , p!9| . (iv) In addition, transport 
equations are usually made tractable by introducing a derivative expansion. We take 
the classical limit of the transport equation, classifying all the terms in powers of h. 
Note that this only partly corresponds to the derivative expansion that is usually 
made. Thus there are two relevant expansions here, which are to be overlaid. We 
treat this double expansion with great care, and explicitly display all powers of h 
explicitly. It is essential at each stage of the expansion, that in addition to the usual 
constraint of energy conservation, chiral symmetry is also preserved. To this end, 
(i) the 1/Nc expansion must be consistently applied and (ii) it is important to see 
that the semiclassical expansion does not destroy the symmetry properties. This 
has been rigorously proven for the lowest order terms in l/Nf.. 

To start with, we examine the terms of 0{{l/Ncf) in the expansion. This 
corresponds to taking the Hartree term only. Such a calculation is usually performed 
self-consistently and this is our approach here also. We then derive the transport and 
constraint equations to 0{{1/ N^)^), and this is done semi-classically, where only the 
leading terms according to powers of h are retained. Note that the aspect of chiral 
symmetry enters into these equations at this point in the property of the subsidiary 
equation, the gap equation, which must be solved concomitantly with the transport 
and constraint equations. 

We next move to terms of order 0{{1/N^^) that give rise to collisions. The 
lowest of these terms contains two interaction lines. According to our philosophy 
of working in a expansion, the self-consistently evaluated Hartree Green func- 
tions must now be employed in this evaluation. In this approximation, we again 
derive the transport and constraint equations, and again use % ordering to decide 
to what level terms in these equations should be accepted. In doing so, we are able 
to demonstrate explicitly, in a fully relativistic fashion, that the Boltzmann-like col- 
lision term is recovered, containing as heuristically expected, matrix elements for 
the elastic scattering processes, qq' qq' and qq' — > qq' , here of course only in 
the Born approximation. In an exact sense, we have rewritten the collision term in 
terms of scattering cross-sections with well defined Feynman graphs, and we have 
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illustrated explicitly which self-energy graphs are required to obtain our result. In 
doing this, we have generalized the non-relativistic treatment of Kadanoff and Baym 
Several other authors have examined the collision term, see for example Refs. 
P,p^. To our knowledge, however, this is the first paper that links the cross-sections 
that occur in the Boltzmann-like equation with Feynman graphs for the scattering 
processes. 

There is, however, a major difference that arises here. By direct calculation, one 
sees explicitly that, even at this level, terms which describe the three body creation 
and annihilation processes, such as g ^ (?^)? as well as vacuum fluctuations, are 
present, but are suppressed by the on-shell requirement of the quasiparticle approx- 
imation. Relaxing this condition would lead automatically to particle production / 
destruction. 

The above calculations have been carried out, in this first exploratory article, 
by considering only one term in the expansion. Clearly the full set is required 
that would completely include the mesonic degrees of freedom. This generalization 
to include the complete mesonic sector is however in itself a complex problem, which 
builds on the concepts and ideas presented here. Again, in this context, care must be 
taken in the h expansion, to see that the properties of chiral symmetry remain intact. 
Using these methods, it should become possible to construct a consistent theory of 
quarks and mesons in which the collision term contains inter alia the hadronization 
processes qq' — > MM' in addition to the usual elastic scattering qq' qq' and 
qq' qq' in the quark sector. The details of the mesonic sector, however, are left 
to a future publication . 

We conclude this section by mentioning some other papers that are relevant, but 
we do not attempt to provide an exhaustive list. Relativistic kinetic equations for a 
Yukawa model were first derived in Ref. Here the difficulties in expressing the 
collision integral in terms of actual Feynman graphs is apparent. A seminal paper 
discussing the Walecka model has followed this see also P]. A transport theory 
for the NJL model, derived from a path integral approach with a semi-bosonized 
Lagrangian was first undertaken by [^. We differ from this in two respects, viz 
(i) our approach is canonical, and treats the theory from a fermionic point of view, 
and (ii) that the terms are ordered according to the double expansion with regard 
to 1/N(. as well as the power counting in h. 
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This paper is organized as follows: in Section II, we review the general non- 
equilibrium transport theory for fermions, derived from a canonical point of view. 
We construct the general transport and constraint equations. In Section III, we 
discuss the quasiparticle approximation, introduce our Lagrange density and explain 
the strategy of the and h expansions. In Section IV, we calculate the mean field 
self-energy and transport and constraint equations in the Hartree approximation. 
We devote Section V to the evaluation of the collision term and demonstrate that the 
classical limit leads us to the Boltzmann equation. Off-shell terms that are dropped 
are seen explicitly. The constraint equation, to this level of approximation, is also 
given. In this section, we present the major new work of this paper. We conclude in 
Section VI. We have, for completeness, included several appendices that give closer 
detail of some of the technical and / or repetitive aspects of the derivations. 



II. GENERAL TRANSPORT THEORY FOR FERMIONS 

A. Equations of motion 

Non-equilibrium phenomena are completely described via the Schwinger - 
Keldysh formalism for Green functions Many conventions exist in the liter- 
ature. In this section, we give ours. For our purposes, it is simplest to use the 
convention of Landau'' 0. In this, the designations + and — are attributed to 
the closed time path that is shown in Fig.l, and the fermionic Green functions are 
defined as 

ihS^x^y) = iTi){x)i}{y)) = ihS''{x,y) 
inS\x,y) = (T'4){x)4){y)) = ihS++{x,y) 
ihS^{x,y) = i^'(p{x)ip{y)^ = ihS~^~ {x,y) 



^ Note that this convention differs from that of Ref. in that those authors attribute 
the opposite signs to the contour branches. There is also another difference that arises in 
that the off diagonal self-energies in this convention of differ from that of Q and by 
a minus sign. 
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ihS<{x,y) = - {i^iymx)) = ihS-^{x,y), (2.1) 

where T is the usual time ordering operator and T the antz-time ordering operator. 
These can be summarized into the compact matrix notation 

S=\ . 2.2 

~ \ s+- S++ ) 

The matrix of irreducible proper self-energies, 

E = , 2.3 

is introduced via the Dyson equation, 

S_{x, y) = S^{x, y) + d'^zd'^wS^{x, w)'E{w, z)S_{z, y) 



= S!^{x, y)+ J d^zd^wS{x, w)Uw, z)S!^{z, y) (2.4) 

where 5° is the matrix of propagators for the non-interacting system. In a standard 
fashion, one constructs equations of motion for the matrix of Green functions. These 
take the form 

{ih - mo) 5(x, y) = a^(5'^(x - v) + J d'^za^T,{x, z)S{z, y), (2.5) 

where 

has been introduced. Note that the entries in the matrices of Green functions, 
Eq.( p.2|) and self-energies Eq. (|2.3|) are amongst one another not independent. They 
are linearly related in a way that is obvious from their definitions, i.e. 

Sapi^, y) + st^i.^^ y) = Sa^i^, y) + s^pi.^^ y)- (2-7) 

It is also useful to introduce the retarded and advanced Green functions, and 
respectively, that are defined as 

inS%{x,y) = 9{xo~y,){{^e.{x),iip{y)}) (2.8) 
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and 

inS%{x,y) = -e{y, - x^){[ilj^{x),i}p{y)]). (2.9) 

Important relationships that connect the Green functions with each other, as well 
as those pertaining to the self-energies, are listed in Appendix A. In particular, it is 
also useful to note that 

S^ = S"" + S"+ (2.10) 

and 

S^ = S— + S+- (2.11) 
link the retarded and advanced self-energies to S*-', i, j = ±. 

B. Transport and constraint equations. 

In this section, we give a brief derivation of the transport and constraint equa- 
tions following the lines of Ref. 0, and writing them in a concise form. It will be 
sufficient for us to examine only one component of this equation, the off diagonal 
element , together with the equation for the hermitian conjugate. This simpli- 
fication arises from the quasiparticle approximation (see below) in which all four 
Green functions 5**-' can be expressed in terms of two distribution functions. The 
equations of motion are 

{ih -mo)„^5^ + (x,?/) = j (x,2;)5^ + (z,?/) + S^^(x,2)S'^+(2;,y)} 

(2.12) 

and 

Sap{x^ y) {-i0y - "^o)^^ = y d^z {-^-^-(x, ^)S;^+(^, y) - S-f^ix, ^)S^+(2;, y)} 

(2.13) 



^We will recover those equations of Ref. |8| in the fermionic sector. Here, however, we 
retain all factors of h explicitly. 
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It turns out to be expedient to eliminate 5*++, S , and S in terms of the 
advanced / retarded forms, together with the off-diagonal elements of the Green 
functions and self-energies, as this will lead us directly to a convenient form for the 
transport equation. Employing Eqs.( [A4| ), ( |A5| ) and ( [A 131 ) - ( |A15| ), one can rewrite 
Eqs.(P7[^) and (P7[B| ) asQ 



{ih ^^-mo)^pSp^{x,y) = j d^z{J:J{x,z)S^^ {z,y) - (x, 2)^^ + (2;, 

z)S^:^iz, y) - z)Sl{z, y)} (2.15) 

and 

Sap{x^y) {-i0y - "^0)^^ = j d^z[-S%{x,z)T.-p^{z,y) + S'^ {x, z)^^^{z,y)] , 

(2.16) 



One now moves to relative and centre of mass variables, u = x—y and X = {x+y)/2. 
A Fourier transform with respect to the relative coordinate, or Wigner transform, 

S{X,p) = J d^ue^'P-'^'^S + - I) (2-17) 

is then performed on the above two equations. Simple rules for the calculation with 
the Wigner transform are given in Appendix B. In general, the Wigner transforms 
of Eqs.(|A5D and (|21^ ) are 



S-+(X,p)AS+-(X,p) - S+-(X,j9)AS-+(X,j9) 
+ S^(X,p)A5-+(X,p) - S-+(X,p)AS^(X,p) (2.18) 



and 



^Note that the more usual form of the first equation would be 

iifi -mo)„^5^+(a:,y) = j d^z[^^S-^ - ^-+S^} , (2.14) 

from which one can immediately see that the right hand side is simply the hermitian 
conjugate of the right hand side of Eq.( 2.16| ). The form that we have chosen, however, 
already anticipates the loss and gain terms that one expects of the collision integral. 
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S-+{X,p){-znr{l^ + ^) - mo} = - 5^(X,p)AS-+(X,p) 

+ ^-+(X,p)AS-^(X,p), (2.19) 

respectively, with 

Now, on subtracting and adding the two resulting transformed equations, one con- 
structs the so-called transport and constraint equations: 

^{7^ ^} + [A S-^X,p)] = /_ (2.21) 



and 



respectively. Here 



2 " ' ^X^' 



y 9X^] + { - mo, 5-+} = /+, (2.22) 



= /coll + + /J, (2.23) 

where the collision term is 

/,oii = S-+(X,p)A5+-(X,p) - S+-(X,p)A5-+(X,p) 

_ rgain rloss (O 0A\ 

— -'coll -'coll ' l^-^^J 

and terms containing retarded and advanced components are contained in 

= -S-+(X,p)A5^(X,p) ± 5^(X,p)AS-+(X,p) (2.25) 

and 

= S^(X,p)A5-+(X,p) T 5-+(X,p)AS^(X,p). (2.26) 

Note that the equations ( |2.21| ) and ( |2.22| ) are general equations that are generic 
for any relativistic fermionic theory. These are already well known, and have been 



derived for the Walecka model in Ref. |@] and rederived for QCD in Ref. |]Tl|], for 
example. 

At this point, one has to specify the particular theory that one wishes to investi- 
gate - in our case, this is the simplest version of the NJL model - and to calculate the 
self-energy terms systematically. In particular, the difficulty that we wish to address 
in this paper, is the evaluation of the collision term in a form that is understandable 
on physical grounds. 
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III. APPROXIMATION SCHEMES AND FEYNMAN-LIKE RULES 



A. Quasiparticle Approximation 



Equations ( p.21|) and (|2.22|) are each a system of coupled equations for the 16 
matrix elements of S~p in spinor space. A complete method of solution of these 
equations can be developed on inserting the spinor decomposition 

- thS-+ = F + 275P + YV, + + ^(T^^^S,,, (3.1) 

for S^^ into the above equation. Similar decompositions for all other quantities 
are required, so that an extremely complicated system of coupled equations results. 
While this guarantees a correct solution from a formal point of view, even in the 
coUisionless system, it leads to a set of 16 x 2 equations that must be solved simulta- 



neously ||2^ , |23[| . Only in limited cases where spatial gradients have been neglected, 
can such a set of equations been solved exactly numerically ||T7| . 

A tremendous simplification arises from the quasi-particle approximation. We do 
this in a fashion which, as will become evident in the following section, represents a 
self-consistent solution within a mean field approximation, i.e. we assume that the 
system consists of fermions that, in interacting only with the mean field, acquire 
a mass m{X) that depends on space and time, and which will be dynamically 
generated due to the breaking of the chiral symmetry of the underlying Lagrangian. 
This is as yet undetermined. The dressed fermions (the "constituent quarks") may 
then collide with each other. The index "if" appearing on the Green functions 
below indicates that they are the mean field ( "Hartree" ) functions. 

P ss' 

+ 6{po + E,)J2vA-p)vs{-p)ft{X, -p)], (3.2) 



and 



S^-{X,p) = -2m^[6{po-E,)Y^Us{p)usip)ff{X,p) 

P ss' 



+ S{po + E,)Y.Vs{-v)vs'{-v)fl'\X, -p)], (3.3) 



ss 
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that introduces the quark and antiquark distribution functions /g*'(X, p) and 
f^^'{X,p) respectively. These functions give the number density per degree of free- 
dom in phase space for on the energy shell quarks, = m^(X) + p^. In each 
case, we have introduced the abbreviation /^^'(X, p) = 1 — /^^'(X, p). The m's and 
w's are normalized according to Us{p)us{p) = 2m5ss' and Vs{p)Vs{p) = —2m5ss'- In 
principle, all spin combinations are possible. We shall restrict ourselves, however, 
to the case /^|' = Sss'fq,q, i.e. we do not allow for spin polarizations. Under this 
assumption, one has 

^^+(X,p) = 2m^^[6{po - E,)fq{X,p) - 6{po + E,)f,{X, -p)] (3.4) 

and 

S^-{X,p) = -2m^^[6{po - Ep)f\{X,p) - 6{po + E,)fq{X, -p)] (3.5) 

using Us{p)us{p) =^ + m and J2s '^s{p)^s{p) =^ — rn. Using this convention, the 
time ordered and anti-time ordered free particle Green functions can be written as 

^ + m 



sr{x,p) 



pZ _ _|_ 

+2m^^[S{po - E,)f,{X,p) + 5{po + E,)f^{X, -p)] (3.6) 



and 

j6 + m 



S^+{X,p) 



p2 _ ^2 



+2m^^-—[6{po - Ep)fq{X,p) + 6{po + E,)f,{X, -p)]. (3.7) 



Note that in terms of the decomposition Eq. (|3.1|) , the decomposition for Sj^ above 
implies that only the scalar and vector components, 

1 777 — 

F = -2-^trV = 2vr/i^[5(po - Ep)f,{X,p) - 6{po + Ep)fg{X, -p)] (3.8) 

and 

V = —F (3.9) 

m 

survive. Here the trace is regarded as being over the spinor variable only. 

Since all four Green functions S^^{X,p) depend only on the two distribution 
functions fq{X,p) and fq{X,p), one needs only to consider the equation for S^^, 
c.f. Eq. ( |2.12|) . Thus the 2 x 16 equations reduce to two. 
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B. The Lagrangian 



In order to proceed further, we must specify our Lagrange density. We choose 
the NJL model form as being 



where G is a dimensionful couphng, [G] =MeV~^. Since the Lagrangian generates 
a non-renormahzable theory, Eq.( p.lO[) has to be supplemented by a regularization 
scheme. We shall use a 3 — cutoff A. For simplicity, we have ignored the flavor 
degree of freedom in Eq.( p.l(]| ). The color degree of freedom is to be implicitly un- 
derstood. This minimal form of the NJL model is chirally symmetric for rriQ = 0, 
and exhibits dynamical symmetry breaking at temperature T = 0. An equilibrium 
study of the model indicates that a phase transition to a chirally symmetric state 
occurs at a finite value of the temperature. By formulating the non-equilibrium the- 
ory for this model, we hope to gain insight into the dynamical consequences of chiral 
symmetry breaking and restoration. The basic problems concerning formulation of 
the transport theory can be studied within the framework of this simpler model, 
and the generalizations to include isospin are obvious. 

From a perturbative expansion in the interaction picture, Feynman-like rules can 
be constructed for the component functions of the matrix of Green functions and 
their self-energies. These follow on considering the same skeleton graphs that one 
would derive for the time ordered equilibrium Green function using an extended 
interaction, and decorating the interaction line with the same sign at each vertex. 
To each interaction line with a + or — sign is attributed a function ^r2iG/hr, with 
F = 1 or for the scalar or pseudoscalar vertex, as is required. External legs have 
a sign attributed to them as well, in accordance with the function to be evaluated. 
Internal lines take all possible combinations of the signs. For example, the Dyson 
equation for the full component S~^'^{X,p) is given diagrammatically in Fig. 2. One 
can see that this diagram and the convention given are in accordance with Eq.( |2.4| ). 
Further examples follow in the evaluation of the self-energy in the following sections. 




(3.10) 
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C. The double expansion in powers of l/N^ and Ti 



Before we evaluate the right hand sides of the transport and constraint equations, 
/_ and /+ respectively, we explain our strategy of a double expansion in powers of 
and h. 

In general, an expansion in powers of the coupling constant is inappropriate for a 
strongly interacting theory. In our theory, Eq.( p.lO[ ), the relevant parameter is GA^, 
which has a value close to two, and thus precludes a perturbative expansion. It is 
therefore customary to perform suitable partial summations of relevant diagrams. A 
perturbative expansion in the inverse number of colors, has been proposed by 

Witten in connection with QCD proper, and such an expansion technique has proven 
extremely useful in connection with calculating observables in the NJL model. In 
this scheme, a diagram containing n interaction vertices and m fermion loops is of 
order G^N^ = {Nc)"'~"^, if one sets GNc = 1, using a symbolic notation. Applied 
to the NJL model, the partial summation at the level 0{{1/Nc)^) is equivalent to 
the Hartree approximation (Fig. 3), and leads to the concept of constituent quarks, 
while to order 0((1/A''c)^), one recovers the Random Phase approximation which 
gives rise to mesonic masses and qq and qq cross sections. 

When investigating the transition of a quantum system to its classical limit, 
one often treats ^ as a small quantity (although this is inappropriate because of 
its dimension) and keeps only the leading orders in h or takes the limit h —>■ 0. In 
the transport and constraint equations, Eqs.( p.21 ) and ( p.22 ), the right hand sides 
I- and J+, respectively depend on h in two ways: (i) The derivative operator A, 
Eq.( p.20[ ), contains h to all orders. Its expansion in h generates the various terms 
of the so-called gradient expansion, and (ii) the self-energies and may have 
overall factors in h. While to order 0{{1/Nc)^), = m, and is directly related 
to the mass, which we treat as a quantity that has a direct classical interpretation, 
we find oc hda/dQ, and is only related via a factor h to a quantity with direct 
physical interpretation. Note that while we use the overall factors in or in 

or S ^ in our h classification, we do not expand the quantities m or da/dQ in 
powers of h. We will detail this strategy in the calculation of the mean field term 
and collision integrals. 
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IV. MEAN FIELD APPROXIMATION 



The first term in the expansion of the self-energy in powers oil/N(. is the Hartree 
diagram of Fig. 3. This is evaluated self-consistently. Only diagonal elements exist 
for Sj:^, i.e. T,Jj^{x,y) = ?/) = 0. One finds by direct calculation that 

— iTj]j^{x,y) = —2ihG (tTiS]j~{x,x) + i'y^tiiSjj' {x, x)i'y5^ 6^{x — y) (4.1) 

while 

- iSj+(a;,?/) = 2ihG (tiiS^-^{x,x) + i-f5tiiS^-^{x,x)i-f5) 5^{x - y), (4.2) 

where tr = tr^tr^ = A^^^tr^. From this self-energy, we may derive two things, (i) the 
gap equation and (ii) the transport and constraint equations at this level. 

Firstly, by examining the form of the retarded Green function for a free parti- 
cle, see Appendix A, one may identify the self-consistently determined mass of the 
particle to be m(X) = mo + S§(X, p = 0) = mo + T.Jf^Xjp = 0) to this level of ap- 
proximation. Because of the 6 functions in Eq.( [4.1| ), S^(X) is in fact independent 
of the momentum variable. Note that Tijf and are scalars in spinor space, 
since the second term in both Eqs. (|4.1j ) and ( [4.2|) vanishes, in view of the quasi- 
particle approximation Eq. (|3.4|) . Now, inserting the quasiparticle decomposition for 
Sjj'{X,p), see Eq. (p.6|) , with m = m{X) into Eq. ( |4.1| ) and performing a Wigner 
transform, we arrive at the gap equation for a non-equilibrium medium, 

d^p 1 

which is a non-linear equation for m(X), and one recovers the result of pT|,|25[ 



mix) = mo + AGNMX) J T^seTx)^^ ' ^"^^'^^ ~ /."(^.P)]. (4-3) 



We may now return to the transport and constraint equations and investigate 
their form. The only contribution from the right hand side of the transport and 
constraint equations, Eqs.( |2.21| ) and (|2.22|) that is non- vanishing, stems from the 



advanced self-energy S;^ that occurs in in Eq.(|2.26|). Closer examination of I^, 



shows that 



2 ' " 

while 

rA _ fsrA 



4{9,^iindi:S],^X,p)} (4.4) 



lt = {^fiiX),S^^X,p)}. (4.5) 
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Here we have made use of the fact that is a scalar and is only a function of X, 
and have expanded A in Eg. ( p.20| ) to retain only the first non- vanishing term. One 
thus has 

if) (9 9^+ fx r)\ 

y {7^ QX>^ } + [i*' ^H^(^'P)] = -^^a,E^(X)a; V(X,p) (4.6) 

and 

Y[7^^] + {i'-"^o,5^+} = {S^(X),S^+(X,p)} (4.7) 

for the transport and constraint equations. From these, one constructs equations 
for the distribution functions fq^q{X,p) by inserting Eq. (|3.4| ) into ( |4.6| ) and ( |4.7| ), 
integrating over Pq over an interval A± that contains Ep{X) or —Ep{X) respectively, 
and performing a spinor trace. Then Eq. (|4.6| ) becomes the Vlasov equation for fq^g, 



P>'d,fq,q{X,p) + m{X)dMX)d';fq,q{X,p) = 0, (4.8) 

where p'^ = Ep{X). We have detailed this derivation in Appendix C, since some care 
must be taken in bringing the equation into this form. The constraint equation, 
Eq.( [4.7|) , which is somewhat simpler to treat, is also discussed briefly in Appendix 
C. On inserting S~^{X,p), and repeating the procedure described above, one finds 

ip'-m'iX))fq,q = 0. (4.9) 

It is interesting to note that, in moving from Eqs.( [4.6D and ( [4. 7] ) to Eqs.( |4.8| ) and 
( [4.9|) , the commutator terms in Eqs.( |4.6|) and (|4.7| ) do not contribute. One sees 



therefore explicitly that the constraint equation ([4.7|) has no multiplicative factors 
of h, while the transport equation contains a simple power of h throughout. This is 
a particular feature of the covariant approach to quantum transport theories which 
has been known for some time [E^,^]. 



The resulting set of equations (|4.8| ) and ( [4.9|) , as is evident, forms the classical 
limit of the set of transport and constraint equations in the Hartree approximation, 
since m(X) is a classical quantity. We note that Eq. ( [4.9| ) is trivially fulfilled for the 
quasiparticle approximation. 
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V. COLLISIONAL SELF-ENERGY 



In this section, we come to the main new aspects of the paper. Here our task is 
to calculate the collisional self-energy for the model, and to manipulate the result- 
ing equations in the semiclassical limit into the form of a Boltzmann-like equation, 
thereby explicitly exposing the weaknesses of the quasiparticle approximation and 
explicitly showing the connection between the self-energies evaluated in this ap- 
proach and the cross-section that appears in the transport equation. In particular, 
we wish to identify the Feynman graphs that are present in the cross-section, and 
their relation to the graphs evaluated for the self-energy. Note that some calcula- 



tions that evaluate the collision term already exist in the literature P-p!0[] . However, 
in some cases, in which the standard spinor decomposition is used [0], the connec- 
tion to cross-sections is not evident and the resulting expressions are unwieldy and 
model-specific. In others, see for example ||^, although it is somewhat simpler to 
connect the result to cross-sections, this has not been done in as general a fashion 
as will be presented here. 

To be specific, we investigate self-energy graphs that contain two interaction 
lines, and which form the contributions to lowest order to the off- diagonal matrix 
elements of E, and which are illustrated for S in Fig. 4. The diagram in Fig. 4(a), 
which we denote as the "direct" graph, is a term of order 1/Nc. The exchange graph 
of Fig. 4(b) has the order Nevertheless, we retain this graph, noting that 

this correction still renders our calculation correct to 0{1/Nc). As will be seen, 
the exchange graph plays an essential role. It leads directly to the cross terms in 
the scattering amplitudes that we will encounter. We note that in a rigorous 1/Nc 
expansion, the propagators occurring in Fig. 4 are now not to be recalculated self- 
consistently, but are to be replaced by the self-consistent Hartree propagators that 
were already determined in the previous section. In addition, the simple use of the 
bare interaction, G, in Fig. 4 at this level is too naive: there are further terms that 
are also of the same order in 1/Nc that must also be included for completeness [|TB,r^ 



Nevertheless, this simplification, which corresponds to a kind of Born approximation 
is a first step in understanding the role of these diagrams in a clear fashion, and 
which paves the way for its generalization to include mesonic intermediate states 
fully [^. The particular highlight of this study is that illustrates explicitly how 
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the relativistically formulated transport theory in the semiclassical hmit results in 
a Boltzmann-like collision term, that can be expressed in terms of a scattering 
cross section which in turn can be evaluated from a well-defined set of Feynman 
graphs. We do this without introducing the added complexity of the complete 
set oil/Nc terms that correspond to the inclusion of mesons. We have nevertheless 
alluded to this extension by our use of notation, labelling the scalar and pseudoscalar 
interaction vertices with a o" or vr respectively. The full extension of this method 
to include the mesonic degrees of freedom will be handled in a separate publication 

We will now evaluate the first term of the collision integral Eq. (|2.24| ) within the 
derivative expansion with A = 1, in order to establish the order in h of this term. 
A further analysis of the remaining terms that occur in J_ can then be made. 

In order to evaluate the collision integral, Eq.( |2.24D , we require the off-diagonal 
self-energies. Consider say . This can be decomposed into three parts, 

= + + S+r^ed(^,P), (5.1) 

and similarly for Here, and refer to Fig. 4 in which all vertices 

have a factor unity or ^75 in spinor space attributed to them respectively. SmJ^ed 
contains the vertices 1 and ^75 in a mixed fashion and is shown graphically in Fig. 5. 
Note that the direct graphs do not appear here, as they are identically zero in the 
quasiparticle approximation. Direct translation of the diagrams of Fig. 4 and 5 in 
x-space can be performed, and then Wigner transformed. This results in 

.2t,2 f d% d% ^^_^^4) 

{2nh)^ (277^)4 {27rhy- 

X [S+-(X,pi)tr {s-+{X,p2)S^'iX,p,) 

-S+-{X,p,)S-+{X,p,)S+-iX,p,)l (5.2) 



K-ix.p) ^ -,G^n^ I - p. + p. - P3) 



while 

X [755+-(x,pi)75tr {S-+{X,p,)j,S+-{X,psh,) 

-755+- (X, pi)755-+(X, P2HS+- (X, p3)75] (5.3) 



and 
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X [755+-(X,pi)5"+(^,P2)755+-(X,P3) 

+S+- {X, pi)755-+(X, p2)S+- {X, p3)75] . (5.4) 

The corresponding expressions for S ^ are obtained from these by exchanging + 
and — . 

Let us now construct the transport equation for fg{X,p). As was indicated in 
Section [TVl , this follows on integrating the transport equation for the Green function 
over the energy variable po on an interval A_|_ containing Ep and taking the trace. 
Applying this procedure to Jcou, we define 

<^coii = tr / dpol^ow 

= J^JPotT {j:^+{X,p)S+~{X,p))- |^^rfpotr(s+-(X,p)5-+(X,p)) 

_ Tg,gain _ Tg,loss (r\ r\\ 

— ^coll ^coU W-^J 



having set A = 1. Now, on taking cognisance of Eq. it is clear that each of 
these terms itself is to be constructed of three components: 

T-g, gain/loss j-g, gain/loss y-g, gain/loss j-g, gain/loss ft^ fi\ 

■^coU — '^coll,cr + "^coU.Tr + '^coU.mixed 

in which the corresponding subindex is placed on E""*" or E"* . In the following 
section, we shall construct the loss term, and do this in detail for the scalar-scalar 
interaction, Jcon a particular. The calculations for the remaining terms then follow 
by analogy, and these are then sketched briefly. 

A. The loss term 

a. The scalar-scalar interation. The scalar-scalar interaction contributes to the 
loss term via 



y-g.loss _ — ^.^ 



27ri 



K {X,Po = Ep,p)J2us{p)us{p)fy{X,p) 



(5.7) 



on inserting Eq.( |3.2| ) into the definition Eq.( |5.5| ), and assuming no spin dependence. 
One may now insert E^~(X, po = Ep,p) from Eq. (|5.2|) into Eq.( p.7|) and multiply 
out the resulting expression. This rather lengthy result contains eight terms. 
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Tg,loss 
"'coll.o- 



ni 4iG^ 



J (27r, 



d'^pi d'^p2 d'^ps 



-(27:h)%E,-pl+pl-pl) 



with 



xS{p-pi+p2 -p3){2nhf 



2Ep^ 2Ep^ 2Ep^ .^^ 



T, = 6{p',-EpJ6 

—ti{uu)p{uu 
T2 = - Ep,)S 

—tii:{uu)p{uu 
Ts = - Ep,)S 

—tr{uu)p{uu 
T^ = 8{p\-Ep,)5 

—iY{uu)p{uu 
T,^5{p'i + Ep,)5 

—tr{uu)p{vv 
T, = 5{p\ + Ep,)5 

—ti {uu)p{vv 
TT^S{p[ + EpJS 

—tT{uu)p{vV 

Ts^5{p\ + Ep,)5 

—tY{uu)p{vV 



P2 - Ep^)^{pl - Ep,^[ii{uu)p{uu)iii{uu)2{uu)^ 

1 (mm) 2 (mm) 3] /g (Pl ) /g (P2 ) (PS ) fq (p) , 
P2 - Ep^)Kpl + -Ep3)[tr(MM)p(MM)itr(MM)2(MM)_3 
(mm) 2 (^;^;) -3] /g (Pl ) /q (P2) /q- ( -Pz)fq (p) , 

pi + Ep^)5{pl - Ep^)[tY{uu)p{uu)itr{vv)-2{uu% 

1 (mm) _2 (mm) 3] (pi ) /g ( -P2 ) /g (P3 ) /g (P) , 

pI + Ep,^)5{pl + £'p3)[tr(MM)p(MM)itr(MM)_2(MM)_3 

1 (mm) _2 (mm) -3] /g (Pl ) /g ( -P2 ) /g- ( -ps) /g (P) , 

pi - Ep^)5ip)l - Ep^)[tY{uu)p{vv)-iiT{uu)2{uu)2, 

_1 (mm) 2{uu)^\ /g ( -Pi ) fq (P2 ) fq (P3 ) /g (p) , 

P2 - -£^P2)^(P3 + -Ep3)[tr(MM)p(MM)_itr(MM)2(MM)_3 
_ 1 (m?7,) 2 (mm) _3] fq ( -pi ) fq {p2 ) /,7 ( -P3 ) f, (P) ■ 

P2 + Ep2)5{pl - Ep^)[tr{uu)p{vv)^itr{vv)^2{uu)3 

_1 (mm) _2 (^i^i) 3] /g-( -Pl ) /g ( -P2 ) /g (P3) /g (P) , 

P2 + -^P2)'^(P3 + Ep^)[tT{uu)p{vv)-itT{vv)-2{vv)-3 
_l(MM)_2(MM)_3]/g(-pi)/g(-p2)/g(-p3)/g(p)], 



(5.8) 



(5.9) 

(5.10) 
(5.11) 
(5.12) 
(5.13) 
(5.14) 
(5.15) 
(5.16) 



where spin indices have been suppressed. Here, the indices "i" mean pj, i — 1,3, 
while the notation "— i" refer to the negative four-momenta, — pj. At first sight, 
this expression appears formidable. It provides us, however, with some physical 
insight. Each term contains a product of a combination of four quark and antiquark 
distribution functions, with fq{p) being present in all of them, being the external 
function under study. One may attribute a diagram to each of these processes in a 
loose sense, by assigning incoming quark/antiquark lines to fq^q and outgoing ones 
to fq^q. In doing so, it appears that quark-quark scattering terms arise from Ti, 
and that quark-antiquark scattering graphs arise from T4 and T7. The remaining 
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graphs correspond to the processes shown in Fig. 6, and which involve particle and 
antiparticle creation / annihilation. Note that time runs from left to right in these 
diagrams. We emphasise that the graphs depicted in this figure are not Feynman 
graphs. The identification of Feynman graphs requires a detailed analysis of the 
coefficient factors, such as that which will be done for Ti, T4 and Tj. In Fig. 6, 
graph (a) corresponds to the process of T3, while (b) corresponds to those of T2 
and T5. The expression for Tg is given by (c), while the vacuum fluctuations of 
(d) correspond to Tq. Note that these all are categorized as off-shell processes: on 
performing the energy integrations over dpi, and dp^, it becomes manifest that 
these terms vanish due to to energy-momentum conservation. Thus, in contrast 
to the non-relativistic case, one notes that the quasiparticle assumption serves to 
suppress particle / antiparticle creation / annihilation, a restriction that may be too 
severe for application to a fundamental theory of quarks and gluons. 

It now remains for us to validate that Ti, T4 and T7 indeed correspond to the 
processes of quark-quark and quark-antiquark scattering and to identify the scat- 
tering graphs in a strict sense. Here we wish to identify actual Feynman diagrams, 
and not simply associate a heuristic diagram to the processes occurring in these 
terms. To this end, we require a detailed analysis of the coefficients. We thus pro- 
ceed by evaluating the energy integrals, and make the substitution pi — > —pi for the 
antiquark state, to find 

rt,ioss ^ AGH^ r d^pi d^p2 d^ps 4 

Ep n J {27inf2Ep-^ {2Tinf2Ep,^ {2Tihf2Ep^ ^ ^ ' 

X [6^{p -P1+P2- Ps) {tT{uu)p{uu)itl{uu)2{uu)3 

-tl{uu)p{uu)i{uu)2{uU3))fq{pi)fq{p2)fq{P3)fq{p) 

+ S^ip-Pl - P2 + P3){ir{uu)p{uu)itr{vv)2{vv)3 

-tr(uu)p(uu)i(vv)2(vv)3)fq(pi)fq(p2)fq(p3)fq(p) 

+ S^{p + P1-P2- P3)itT{uu)p{vv)itr{vv)2iuu)s 

-tr {uu)p{vv)i{vv) 2{UU) 3) fq{pi)fq{p2)f q{p3)fq{p)]- 

(5.17) 

This can be simplified slightly, if one makes the substitution pi ^ p^ in the last 
term, enabling one to combine the second and third terms. At the same time, one 
notes that there is a symnmetry in the first term in the variables pi and p^. We 
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employ this to rewrite the first term also. One then has the form 
rt,loss 46*2^3 f d^p^ d^p^ d^p3 



jq,loss / ('9 ^i^l^ 

~ ~Yp~h~ J {2nny^2Ep^ {27cn)^2Ep^ {2nhf2Ep^ ^ ^ ' 
X +P2 -p3)^[tr(nn)p(MM)itr(nn)2(uM)3 

+tr ('U'u)p (mm) str (mm) 2 (mm)! — tr(MM)p(MM)i(MM)2(MM)3 

-tr(MM)p(MM)3(MM)2(MM)i]/g(pi)/g(p2)/g(P3)/<7(p) 
+ -P2 +P3)[tr(MM)p(MM)itr(MM)2(MM)3 

+tr(MM)p(MM)3tr(MM)2(MM)i - tr(MM)p(MM)i(MM)2(MM)3 
-tr(MM)p(MM)3(MM)2(MM)i]/g(pi)/5(p2)/g(P3)/9(p)} 

(5.18) 

Now, by directly examining the quark-quark scattering processes, see Fig. 7, one sees 
that the first two terms that arose from the direct graph of Fig. 4(a) correspond to the 
t and M channel exchange diagrams respectively, while the remaining terms that arose 
from the exchange graph of Fig. 4(b), are constructed from the product combinations 
of the t and u channel matrix elements. The individual matrix elements for these 
processes are 

- (gg ^ qq) = '^uy^ulul (5.19) 
-iM^iqq ^ qq) = ^u'^u^ulu^^. (5.20) 

Similarly, in the quark antiquark channel shown in Fig. 8, the s and t channel graphs 
are represented together with their products. Their individual matrix elements are 

2tG 

T 

-zM*(gg ^ qq) = ^uy;v§vl (5.22) 

Squaring these matrix elements and comparing the result with the terms in 
Eq.( |5.18D , one arrives at the simple form 

jloss,a ^ _^ f d^Pl d^P2 d^P-i .g^x4 

« Ep J {2Tih)nEp^ (27r^)32Ep, {2-nh)nEp.^ ^ ' 

X {6\p - Pi + P2 - P3)\Y: \M^q^Jp2 l3)\'fMfMMp3)fqip) 
+ 5\p -PI-P2+ Ps) E IK^^.,-(P3 - 12)\'fMfMf,ip;)f,ip)}, 

(5.23) 



- ^M^(gg ^ qq) = n^^"^^ < (5-21) 



23 



where |M^^__^gg| = |M*(gg — > qq) — M^{qq qq)\ and IM^^^^^I = \M^{qq — > qq) — 
M^{qq —>■ qq)\ contains the relevant s, t and u channels, and the sum refers to all 
spin variables. 

In this calculation, one has explicitly observed that in the quark-quark scattering, 
the direct graph of Fig. 4(a) leads to the squares of the matrix elements of both the 
t and u channels, while cross terms are attributed to the exchange graph. Similarly, 
for the quark- ant iquark processes, the direct graph of Fig. 4(a) leads to the squares 
of the matrix elements for the s and t channels, again with the cross terms arising 
from Fig. 4(b). One thus sees the necessity of retaining both of these graphs in order 
to derive the sum of different channels in the cross section. 

b. The pseudoscalar-pseudo scalar interaction. The interaction graphs of Fig. 4 
that contain only the pseudoscalar ^75 at each vertex contribute to the collision term 
via 



y(J,l0SS _ ^.^ 



27ri 



'coILtt 



2Ep 



(5.24) 



with E~"*" given in Eq. ( p.3| ). This term may be dealt with in the same fashion as the 
scalar term. Explicit details are left to the Appendix D. After some manipulation, 
one finds 

Tioss,^ f d^Pi d^P2 d^P3 ^o^^^4 



Ep J {2TTh)^2Ep^ {2TTh)^2Ep^ {2'Kh)'^2Ep,^ 

X {5\p -p,+p2-p,)\y: iK™(p2 - m'fMfMiMUp) 

+ 5\p -P,-P2+ Ps) E IK™(P3 - 12) r/,(Pi)/^(P2)/,(P3)/,(p)}, 



(5.25) 

where {M^q^^^l = \Ml{qq qq) - M^{qq qq)\ and W^g^^^l = \M^{qq ^ qq) - 
Ml[qq — »• qq)\ contain the relevant s, t and u channel scattering processes that are 
given in Figs. 7 and 8 with the vertices ^75 everywhere. 

One can now combine Eqs.( ^.23D and ( ^.25|) . These are added. In addition, it 



is useful to relabel the variables 2 3 in the quark- ant iquark scattering terms, so 
that one may extract a common 5-function. One then has the combined terms 



tIoss,^ jq,\oss . jq, loss 

•^coll,CT+7r ~ ■^coU.o- + ■-'con,n 



J i2nh)^2Ep, i2nn)^2Ep, {2.h)nEp, ^'"^^ ^ (^-^i+^^-Ps) 



Ep J {2Tin)^2Ep^ {2TTh)^2Ep^ {2nn)^2E. 
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X {^(E IK,->,.(P2 - 13)P + E |M;,^^^(p2 ^ 13)n 

^fgiPl)fqiP2)fgiP3)fqip) 

+ (E IK™(p2 - i3)p + E |m;,^.,-(p2 ^ i3)n 

xfgiPl)fgiP3)fgiP2)fgip)}. (5.26) 

From this form, it is evident that both scalar and pseudoscalar interaction terms 
combine to give the individual matrix elements squared. What is still missing, 
however, is the mixed self-energy. This is examined briefly in the following section. 

c. Mixed interaction. The mixed self-energy, given in Eq.( ^.4|) contributes to 
the loss term via 



T5,ioss ^ 27ri 

coll, mixed 2,E 



Siixed(^.PO = Ep,p)Y,Us{p)Us{p)fg{X,p) 



(5.27) 



Dealing with this term as with the scalar interaction, one finds a series of terms, 
the details of which are to be found in Appendix E. One can identify these in terms 
of the amplitudes of the scalar and pseudoscalar scattering matrix elements, see 
Eqs. dCTp , {^^), ([02]) and (0) and one finds 

rt.loss _ f d^Pl d,^P2 d^PS /o ^x4j 



+ E M^M** + E M* Mr + E Km:* + E km! * 



= J {2.hY2E,, {2.Kf2E,^ {2.nf2Ej ^''^~^ -Pi+P^- Ps) 

X {1[E M^Mi* + E M*Mr + E MiMr + e m-m**],,^,, 

x/g(Pi)/g(P2)/g(P3)/g(p) 

^fg{Pl)fg{P3)fqiP2)fliP)}- 

(5.28) 

Thus, the final result for the collision term, J^on^^, constructed by adding Eq. ( p.28| ) 
to Eq.dOBI), is 

^-^^ = -E-J {2.hf2E,^ {2.hf2E,^ i2nh)^2Ej^''^^ ^ (p - Pi + P2 - Ps) 

X{^E I^..-..(P2 - m'fgipi)fgip2)fgiP3)fgip) 

+ Y.\Mgg^gg{p2 ^ IS)]' f,{pi) fM fM f,iP)} ^ (5-29) 

where the matrix elements are now complete in both the quark-quark and quark- 
antiquark sectors, Mg^^qq = M^^^^^ + M^q_^^g and Mqq^qq = M^^^^g + M^g^^^. We 
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comment once again that although this resuh has been derived using a constant 
interaction or in Born approximation, it will also hold in an extended form, where 
the interaction proceeds rather via a momentum dependent coherent mesonic state 

m. 



B. The gain term 



Construction of the gain term proceeds exactly as in the loss case with the 
evaluation of 



Tq.gam 



2Tci 
2E„ 



-tr 



S +(X,]9o = Ep,p}J2usip)usip)fqiX,p} 



(5.30) 



for the three different forms of S We do not repeat this procedure here, but 
simply quote the final result. One has 

Tih r d^pi d^p2 d^P3 



x(j,gam 
>^coll 



Ep J {2'Khf2Ep^ {2'Khf2Ep^ {2T:hf2Ep^ 

+ J2\M,,^,,ip2 ^ m'fMMps)Mp2)Mp)} 



i2Ti?iY5\p-p^+P2-p3) 



(5.31) 



with the complete quark-quark and quark- ant iquark scattering matrix elements as 
defined in Eq.(^:29D. 



C. The transport equation 

In this section, our aim is to construct the transport equation that now includes 
the self-energies represented in Figs. 4 and 5. In order to do so, we first combine 
the results of Section |V A| and |VB| to construct the collision integral J^^u that was 



defined by Eq.( |5.5| ). This energy integrated collision term can then be written as 



eJ ^ J {27ihf2Ep, 



Jlii = -— IdQ 

+ E I^m(p2 - m\Upi)Up2)h{p^)Up) - Upi)Up,)h{P2)Up)}^ 

(5.32) 
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where 



rlQ = i2.nmp + p., - - p3) p^ll^ (5.33) 

is the Lorentz invariant phase space factor. One can, however, reexpress the right 
hand side of this equation in terms of the differential scattering cross section. One 
can see this simply if one considers centre of mass variables: noting that 

dQ = ^ . l^^l dn (5.34) 

and 

da 1 \pi\ 



dVL GAtt^'s \p: 



\M\\ (5.35) 



'2 



where \M\ is the spin and color averaged scattering matrix element squared. One 
thus has 



JLi = §r- dn ^^-f^\v, - V2\2E,2E, 



P2 



Ep J J {2Ti?if2Ep^ 

X {i^|^,_^^(p2 ^ imMlMfMUp) - Kpi)fMfMUp)] 



(5.36) 



for the collision term. We have thus obtained the final expression for the collision 
integral that looks formally similar to that expected in the semi-classical limit. In 
addition, we note that this term is proportional to h, explicitly demonstrating our 
a priori claim of Sec. Ill C that terms containing off-diagonal matrix elements of S 
multiplied with 5* are proportional to h. 

In view of this, we may examine the right hand side of the transport equation 
Eq.( p.2lD that contains /_ of Eq. ( 12!2^ ) . We require the energy integrated and traced 



form of J_, i.e. 

Jl = Jcou + J'-"" + J'-'', (5.37) 

with 

j,_,A/R dpol^'"". (5.38) 
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We are required to consider J^'"^ and J'^^ (or alternatively and ) also only to 
this first order oi h for this term in the self-energy. Examining Eqs. (|2.25| ) and ( p.26| ), 
one notes that J'^^ and J^!^ then vanish identically, due to the cyclic property of 
the trace. Higher orders would be of order 0{h^), and may therefore be discarded. 
Thus the transport equation Eq. (|4.6|) now becomes modified to read 



if) r)'^~^(X r)\ 

y {7^ -^y^} + [i*' 5-+(X,p)] = -ind,^UX)d^,S-''{X,p) + /,ou, (5.39) 

with Jcoii related to Jcoii via Eq. ( |5.38D and Jcoii as given in Eq.( p.3(j| ). It becomes 
clear, on integrating over the positive energies po and performing the spinor trace. 



^/ dp, tr{7^ ^^7lf'^^ > = -^^/ dp,d,^jj{X)dt;S-+{X,p) + J,,,, (5.40) 
2 JA+ aA^* JA+ 



that all terms that survive this procedure, are directly proportional to and that 
the resulting transport equation 



p^d.f, ( X,p)+m{X)d,m{X)d^J,{X,p)=Kj dn I ^-f^^\v^-v,\2Ep2Ep, 

+^Im(p2 - imMfMfMm - iMfMfMMp))}, 

(5.41) 

is the classical result. Here this equation has been cast into a form that is reminiscent 
of the non-relativistic Boltzmann equation To make this analogy exact, one 

should regard this equation rather in terms of the quark density distribution, that 
is given as nq{p,X) = 2Ncfq{p, X). It is thus apparent that the lowest order term 
in the derivative expansion of Jcou gives the relativistic analogue of the classically 



expected collision term [^. By contrast, the first order term was required in the 
local self-energy in deriving the force term from the Hartree approximation. Finally, 
we comment that the transport equation for the antiquark distribution function is 
simply obtained from the aforegoing by replacing g g everywhere. 

Before concluding this section, we make a comment with regard to the relevance 
of chiral symmetry to Eq. ( |5.4lD . Formally, any arbitrary interaction will lead to a 
similar form as has been given in Eq. ( p.41[ ) for the transport equation, and therein 
lies the generality of this calculation. However, a different self-energy would de- 
scribe the mean field part or Vlasov part. For the NJL model, this amounts to 
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a dynamically generated quark mass m(X) appearing, that arises because of the 
symmetry breaking. Note that in the Wigner-Weyl realization of the symmetry, 
this term would vanish identically. In addition, the issue of chiral symmetry reflects 
itself again in the cross-sections that are to be calculated for Eq.( |5.41| ), see Ref. |p6|| . 



D. The constraint equation 

The constraint equation, Eq. (|2.22| ), can also be examined in light of the h expan- 
sion. One notes that it is intrinsically of one order less than the transport equation 
in terms of an h expansion: the semiclassical form of Eq. ( p.22| ) arises from struc- 
tures. Thus, in extending the Hartree equation to include coUisional self-energies, 
one may immediately disregard the coUisional term Jcou, as well as since they are 
one order higher in h than the Hartree approximation. An additional correction to 
the mass could at most arise from but which represents 0{l/Nc) corrections to 
m{X) = mi{{X) calculated via Eq.( ^.3| ). In the absence of this term, the constraint 
equation retains the form already given in Eq.( |4.y| ) 



VI. SUMMARY AND CONCLUSIONS 



In this paper, we have presented a general derivation for transport and constraint 
equations for a fermionic system from a canonical point of view along the lines of Ref. 
p|. We have evaluated these for the single flavor NJL model. Here two expansion 
schemes are overlaid. The strong coupling theory is dealt with via an expansion 
in the inverse number of colors, l/N^, in what is by now standard practice. We 
make an expansion in h, whch only partly is identical with the gradient expansion. 
We have evaluated the lowest order self-energy term in the 1/Nc expansion, and 
constructed transport (Vlasov) and constraint (mass-shell) equations. In addition, 
that one is working with a chiral model evidences itself in that the gap equation must 
be solved in addition to the transport and constraint equation emerges, and must 
be solved together with these. It is demonstrated explicitly that differing orders in 
h are required for these two equations 0, the constraint equation being of order 



^Although this feature is known in the Uterature [22|, it is not often made expUcit [p4| 
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while the transport equation has order h^. 

We then examine a term of the next order in the l/N^ expansion. We treat the 
simplest one that gives rise to a non-vanishing collision integral, and show explicitly 
that it is of order h, and therefore is necessary for constructing a consistent transport 
equation. By the same token, it has no relevance for the constraint equation, for 
which only additional mass changes can occur. 

In the derivation of the relativistic Boltzmann equation, we have gone to some 
trouble to indicate explicitly that the on-mass shell condition restricts the collision 
integral to contain only elastic scattering processes, and suppresses pair creation 
and annihilation amplitudes q ^ {(l<i)<i and also vacuum fluctuations that would 
otherwise be present. The derivation of the collision integral that we have given 
here is general for all relativistic two body interactions - minimal additional ma- 
nipulations are required for a momentum dependent interaction. Once again, the 
addition feature of the chiral symmetry manifests itself in that one has to solve a 
gap equation concomitantly. In addition, the cross-sections are a function of the 
dynamically generated quark mass, and therefore also reflect the chiral symmetry 
aspect of the problem. 

In our calculation, we have also been able to identify the associated t and u 
channel matrix elements required to this level of approximation for quark-quark 
scattering, as well as the s and t channel amplitudes for quark-antiquark scattering, 
i.e. we are able to explicitly attribute actual Feynman graphs, and not simply 
heuristic diagrams, to each of these processes. 

In this paper, we have not made any attempt to include the Fock term, which is 
also of order 0{1/Nc). This is deliberate. Since the next stage in this study involves 
the incorporation of the mesonic sector, of which the Fock term forms a part, it has 
no relevance here at this stage and can at most modify the quasiparticle mass arising 
from the Hartree term. It is important to include all diagrams to 0{1/N^, thereby 
incorporating the mesonic sector completely. This in itself is a complex task, and 



contains several other difficulties. This will be the subject of a future paper pO 
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APPENDIX A: GREEN FUNCTIONS 



In this section, we list relations among the various Green functions and the 
self-energies. They can, in almost all cases, be simply derived from their definitions. 

The causal and anticausal Green functions, see Eq. (|2.1| ), are related to S^'^ and 
S^^ via the relations 



and 



(Al) 



y) = ^(^0 - xo)S^i^ (x, y) + ^(xo - yQ)S^^{x, y). 



(A2) 



The Green functions are not independent. They are linearly related in a way that 
is obvious from their definitions. 



Sap {x, y) + 5+3+ (x, y) = S^flix, y) + 5+^ (x, y) 

= Sal3{x,y) = -J^{[i^a{x),i^p{y)]) 



(A3) 



and which defines the Keldysh Green function. On the other hand, the retarded 
and advanced Green functions that were introduced in Eqs.( p.8|) and (2^) can also 
be related to the S"^^ , i,j = ±, via 



Sapix,y) = 6{xo-yo) {S^p {x,y) - S^^{x,y)} 

= S^p{x,y) - S-^{x,y) = S^p{x,y) - S^^{x,y) 



(A4) 



and 
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= S^fsi^, y) - y) = y) - 5+3+ (x, y). (A5) 

In addition, the following relations are evident: 

Sap{.x,y) - S^p{x,y) = S;^p{x,y) - S-^{x,y) (A6) 
SU^.y) + S%{x,y) = S-^-ix,y) - 5++(x,y). (A7) 

It is often useful to alternatively consider the matrix of Green functions 

S'=('''], (A8) 

that can be obtained from by a transformation 

^ = R-^S R , (A9) 

with 

The transformed self-energy matrix S' is 

yK yR 

^'=11. 1 • (All) 







with 





= 


+ s 




= j:— 


+ s 




= j:— 


+ s 



(A12) 
(A13) 
(A14) 

This can be seen directly, on using the equation 

S++ + S— = -(S+- + S-+), (A15) 



which, in itself can be derived via Eq.(|A3|). The Dyson equation for the transformed 
matrix is 
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Sl.{x, y) = 5:^°^ (x, y)+ j d^zd^wS!,^^^ (x, w)EL{w, z)SL{z, y) (A16) 

leading to equations for the retarded and advanced Green functions that are indi- 
vidually closed: 

= 5?;^^°^ (^,2/) + / d'zd'wS''^f'\x,w)T}'^i\w,z)SX{z.y) (A17) 

{ih - mo)^f^Sf:^^{x,y) = Sa-yS'^ix - y) + J d'^zi:^f{x,z)Sf:f"{z,y). (A18) 

The equation of motion for the Keldysh Green function is found on multiplying 
out Eq. ( |A16|) and using the fact that (i — mo)^^*^*^-* = 0. One then has 



{th - mo) S"" = J + (A19) 

i.e. an integro-differential equation, whose interpretation should also lead to the 
generalization of the Boltzmann equation |^. This, together with Eq.( |A18D form a 
complete set of equations. Usually, the solution of either or is required as 
these are hermitian conjugates of one another. 

We now write down the solutions for the free particle case. It is simplest to 
evaluated S^'^. A direct evaluation leads to the form 



S^'^ip) = , ^ + 7- ■ (A20) 



When a self-energy is included, one has 



Using Eqs.O, Q and Q, one may verify that S—{X,p) and S++{X,p) take 
the forms that have been given in Eqs.(p.6|) and ( |3.7| ). 



APPENDIX B: WIGNER TRANSFORMS 

The transformation rules for applying the Wigner transform to the functions 
required in the equations of motion are 



33 



d^fi^,y)^i^^ + ld'x)fiX,p) (Bl) 

d!if{x,y)^H^ + ^d^x)fiX,p) (B2) 

m9{x,y) /(X)exp 9{X,p) (B3) 

mgix,y) /(X)exp (B4) 



These relations can be derived simply from the definition for the Wigner transform, 
Eq.( |2.17[ ) and retaining only first order derivatives in the expansion need be kept. 



APPENDIX C: VLASOV AND CONSTRAINT EQUATIONS 

In this Appendix, we derive the Vlasov equation from the Hartree approximation 
to the self-energy in the transport equation, and the constrained equation in this 
approximation. Starting with the transport equation, Eq.( 



it) (9 9"+ fx n) 

-j{Y,^^^} + [^^,s-Hx,p)] = -^n^,J:\x)^^s-^x,p), (ci) 



we insert Eq.( p.4|) for S ^ for the Green function. The anticommutator on the left 
hand side of Eq. ( |C1|) does not vanish, whereas the commutator does after a spinor 
trace is performed. One finds 
d p^ 



' 2Ep 



[6{po - E,)f,{X,p) - 6{po + E,)fg{X, -p)]) 



+ miX)d,miX)d^i^JKPo - E,)f,{X,p) - 5(po + Ep)f,{X, -p)]) = (C2) 

where a spinor trace has been performed, and we have made use of the fact that 
d^T,'^{X) = df^m{X), since mo is a constant. Now this equation is integrated over 
a positive energy range. In this fashion, only the quark distribution function can 
survive. One thus has 

+ m{X)dMX) I dpod^^S{po - Ep)f,{X,p) = 0. (C3) 
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It is useful to write the Lorentz indices out explicitly. Using the notation i = 1..3, 
this reads 

+^(X)9.m(X)9;(A^|l^) = 0. (C4) 

Note that no term proportional to m{X)dQm{X) occurs, as this term contains a 
vanishing surface integral. If one explicitly performs the differentiations indicated, 
noting that V^^"^ = -{m{X)Vm{X))/E^ and VpE'^ = -p/E^, one finds 

-^f,iX,p) + v- V/,(X,p) - VEp ■ VMX.p) = 0, (C5) 

where we have written v = p/Ep. This is the Vlasov equation, and it is manifestly 
no n-co variant. One can, however, extract a factor Ep^ from the above equation, 
and write 

^(p'^dMX,p) + m{X)dMX)d^f,iX,P)) = 0, (C6) 

giving it a covariant appearance. Of course, in the above expression, p° = Ep, while 
there is no po dependence in the distribution function. It is this equation that has 
appeared as Eq.(|4.^) in the main text. 

The derivation of the constraint equation is far simpler. Consider Eq. ( |4.7] ), and 
insert the quasiparticle form Eq.( |3.4D for S^^. The commutator, which is the first 
term of this equation, gives no contribution upon taking the spinor trace. The 
remaining anticomutator terms can be combined, to give 

tT{^-m{X),S-+} = (C7) 

since m{X) = niQ + T,^{X) = rriQ + in this approximation. Integrating on 

A+ leads directly to the equation 

{p'-m\X))f,{X,p) = (C8) 

that was quoted in Eq.(^.9D. 



APPENDIX D: PSEUDOSCALAR COLLISION TERM 

In this appendix, we provide the detail required to construct the collision term 
from Eq. (|5.24]) with the self-energy as calculated via Fig. 4 using the pseudoscalar 
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vertices Z75. The self-energy has been given in Eq.( ^.3|) . On inserting this into 
Eq.( |5.24D and muhiplying out the result, one again finds inter alia off-shell pro- 



cesses of the form shown in Fig. 6, where the interaction vertex is now pseudoscalar. 
Retaining only those terms that survive energy-momentum conservation, and per- 
forming the same manipulations that led to Eq.( p.l7D , one finds the analogous result 



rt.ioss _ TTi iC^i^ f d^pi d?p2 d^Ps 



~ J {27Thf2Ep^ {2Txhf2Ep^ {2T^Tif2Ep^ ^ ^ ' 

X {S'^iP -P1+P2- Ps) [tr(uu)p75(uu)i75tr(uu)275(^i^i)375 

-t r (mm)p75 (mm) 1 75 (mm) 275 (mMs) 75] /g (Pl ) {P2 ) fq (P3 )fq{p) 
+ 5'^{p-pi -P2 +P3)[tr(MM)p(MM)ltr(MM)2(MM)3 

-tr(MM)p75(MM)i75(MM)275(w)375]/g(Pl)/9(P2)/g-(p3)/g(p) 
+ 6^{p + Pl -P2 -p3)[tr(MM)p75(MM)i75tr(MM)275(MM)375 

-tr(MM)p75(MM)i75(MM)275(^i^i)375]/g(Pl)/9(P2)/g(p3)/g(p)}- 



(Dl) 



Note that the terms with two traces arose from the direct graph Fig. 4(a) while the 
term with one trace originated from the exchange graph Fig. 4(b). now, once again. 
Fig. 7 may be examined with pseudoscalar vertices, leading to the t and s channel 
forms, 

- zMiiqq ^ qq) = '^u^^i^.^.TPu^^u'^^ {i^Y ^' 4 (D2) 
-iM:{qq qq) = '^u^^i^.^.Y^^u^^ {i^^ ^' u^^ (D3) 

for quark-quark scattering. The squares of these terms correspond to the first 
two terms of Eq . ( P 1|) , while the latter two terms are given by the combination 
X](M"M^* + M^M"*). In a similar fashion, the quark- ant iquark channels from 
Fig. 8 are 

- tM^^{qq ^ qq) = '^uUil,)''^v^v^' (t^,)"'^' u^^ (D4) 
-iMliqq ^ qq) = '^u'^(^^,)'^Pu^^v^' (^^,)-'P'v( (D5) 

which combine in the same way to yield the terms that form the sum multiplying 
the second (5-function. One thus has 
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Ep J {27ih)^2Ep^ {2nh)^2Ep^ {2Tihy^2Ep,^^ ' 

+ -p,-p,+ p,) ^ |M;^^^^-(p3 ^ 12) |7,(Pi)/^(P2)/,-(P3)/,b)}, 

(D6) 

where \M^q_qq\ = \Ml{qq qq) - M^{qq qq)\ and \M^g_gg\ = \M^{qq ^ qq) + 
Ml{qq — > qq)\ again contain the relevant s, t and m channels. This was quoted as 
Eq.dH). 



APPENDIX E: MIXED COLLISION TERM 

In this appendix, the contribution of Sj^j^^^ to the collision term is examined. As 
before, and 5*+" are inserted into Smixcd given by Eq. (|5.4|) , and the expression 
for JcoHmixed Eq. (|5.27|) is multiplied out. One may directly identify the energy- 



momentum conserving terms in this case to be 

.g.loss TTi AC^i^ f d^pi d^P2 rf^Ps 



jq,loss " " / fO 1=i\^ 

coii,mixed - y (^27Th)^2Ep^ {2Txhf2Ep^ {27cnf2Ep.^^ ^ ' 

X {S{p-pi +P2 -p3)[tr(MM)p75(MM)i(MM)275(MM)3 

+tr (uu) p (uu) 175 (uu) 2 (mm) 375] {pi ) /s (^2 ) fq (ps) /</ (p) 
+ -p2 +P3)[tr(MM)p75 (mm) 1(^)275(^)3 

+tr(MM)p(MM)i75(w)2(w)375 + tr(MM)p75(w)3(OT)275(MM)l 
+tr{uu)p{vv)3-f5{vv)2{uu)a5]fg{pi)fg{p2)fg{p3)fq{p)}- (El) 

A slight simplification can be introduced on letting 2 3 in the second term. One 
has 



loss 



ni AG'^i^ f d^pi d^p2 d^ps 



coll, mixed 



Ep h J {2nh)^2Ep^ {2TiKf2Ep^ {2rrn)^2Ep^ 
X {2TThY6{p -P1+P2- Ps) 

X [tr(MM)p75(uu)i(MM)275(MM)3 + tr (Mn)p(MM) 175 (uu)2 (mm)375] 
x/g(pi)/2(p2)/g(p3)/<?(p) 

+ [tr(uM)p75(uu)i(w)275(w)3 + tr(MM)p(uu)i75(i;i;)2(w)375 
+ti{uu)p'y5{vv)3{vv)2'y5{uu)i + ti{uu)p{vv)3'y5{vv)2{uu)i'j5] 

Xfq{Pl)fqiP2)fqiP3)fqiP)}- (E2) 
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It is now essential to identify these terms. In order to do so, we note that, for 
the scattering processes of Fig. 7 that include both t and u channels, the complete 
scattering amplitude is 

M„ = M* - M: + M* - M:. (E3) 

Squaring this leads to 

\Mgg^g,\^ = | ^ + | | ^ + M ^ * " 



- M^Ml* + M:^M^*. (E4) 



The first two terms of this expression have already been generated in Eq.( 5.26| ). 
Performing the sum over spinor indices, one sees from Eqs. (|5.19|) , ( ^.2(]| ), (P2|) and 
), that 

Y: MiMl* = = ^ M^Mr , (E5) 



but that 



Y M^Ml* = —ii{uu)p{uu)i-i^{uu)2{uu)3-i^ (E6) 



Now, direct inspection of Eq.(|El|) indicates that all the terms of Eq.(EB) are present, 
if one extracts a factor 1/2 and makes use of the symmetry in the variables 1 and 
3. The qq scattering term can be viewed similarly in its s and t channels for tt and 
cr, and the terms may all be identified. One thus finds the result for JcoU mixed that 
was quoted in Eq.( p.28| ). 
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FIGURE CAPTIONS. 

Fig.l : Closed time path on which the Green functions are defined. 

Fig. 2 : Diagrammatic expansion of the Dyson equation for the Green function 
S++. 

Fig. 3 : Self-consistent Hartree approximation to the self-energy. 

Fig.4 : Direct (a) and exchange (b) graphs that contribute to E""*" and which 
contain two interaction lines. 

Fig. 5 : Mixed graphs that contribute to S^j^g^j and which contain two interaction 
lines. 

Fig. 6 : Processes of the form q and {qq)q — > g as well as vacuum fluctu- 

ations that are suppressed in the quasiparticle approximation. These are heuristic 
graphs and are not Feynman graphs. They should be read from left to right. 

Fig. 7 : The t and u channel Feynman graphs for elastic quark-quark scattering. 

Fig. 8 : The s and t channel Feynman graphs for elastic quark- antiquark scatter- 
ing. 
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